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Abstract 

Helicases, involved in a number of cellular functions, are motors that translocate along single- 
stranded nucleic acid and couple the motion to unwinding double-strands of a duplex nucleic acid. 
The junction between double and single strands creates a barrier to the movement of the heli¬ 
case, which can be manipulated in vitro by applying mechanical forces directly on the nucleic acid 
strands. Single molecule experiments have demonstrated that the unwinding velocities of some he¬ 
licases increase dramatically with increase in the external force, while others show little response. 
In contrast, the unwinding processivity always increases when the force increases. The differing 
responses of the unwinding velocity and processivity to force has lacked explanation. By general¬ 
izing a previous model of processive unwinding by helicases, we provide a unified framework for 
understanding the dependence of velocity and processivity on force and the nucleic acid sequence. 
We predict that the sensitivity of unwinding processivity to external force is a universal feature 
that should be observed in all helicases. Our prediction is illustrated using T7 and NS3 helicases as 
case studies. Interestingly, the increase in unwinding processivity with force depends on whether 
the helicase forces base pair opening by direct interaction or if such a disruption occurs sponta¬ 
neously due to thermal fluctuations. Based on the theoretical results, we propose that proteins 
like single-strand binding proteins associated with helicases in the replisome, may have co-evolved 
with helicases to increase the unwinding processivity even if the velocity remains unaffected. 
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I. INTRODUCTION 


Hclicases are molecular motors that can translocate along single-stranded (ss) nucleic 
acids (NAs) and can also unwind double stranded (ds) NAs when confronted by a single 
strand-double strand junction DM]- These remarkable motors, found in all organisms, are 
involved in a variety of cellular functions ranging from DNA and RNA replication, remod¬ 
eling of chromatin, translation and other aspects of nucleic acid metabolism [SHU], and are 
classified into six families based on their sequence PUGEj. Various ensemble assays have 
been developed P2HZE! to decipher how they step on single strands of nucleic acids, how 
they interact with single strand-double strand junctions, how often they dissociate from 
their track. Single molecule manipulation techniques, such as laser optical tweezers (LOT) 
and Magnetic Tweezers (MT), have also been used to probe the kinetics of stepping and 
nucleic acid unwinding at the level of individual helicase molecules [22h25] . These exper¬ 
iments are particularly suited to obtain quantitatve measurements of helicase processivity 
and unwinding velocity under force. From such measurements, one can obtain insights into 
the mechanisms of helicase function, and how various hclicases differ from each other. 

In a number of single molecule experiments, a variety of interesting observations have 
been made on the response of helicase velocity and processivity to force. Processivity is a 
measure of the average number of base-pairs unwound by the helicase in a single binding 
event (see below for a more detailed discussion of various other definitions). Dessinges et al. 
[21], who used MT to follow the unwinding of dsDNA by the UvrD helicase, showed that the 
unwinding velocity depends only weakly over a broad range (3-35 pN) of applied forces on 
the ds termini of the DNA. However, the unwinding processivity was much larger (265 bp) 
in their experiment compared to previous results (45 bp) from zero-force ensemble assays 
[18] , suggesting that force possibly plays a part in enhancing the processivity of 
UvrD. It should be kept in mind that the stoichiometry of helicase attachment 
to DNA was not explicitly determined in |23j- This is particularly relevant 
because enhanced processivity is achieved by UvrD dimers | ]26L 12T J. In sharp 
contrast, Johnson et al. [28] discovered using LOT experiments, that both the unwinding 
velocity and processivity of T7 helicase are highly tension dependent. Indeed, the unwinding 
rate increased by an order of magnitude when the tensile force applied to destabilize the 
ss-dsDNA junction was increased from 5 to 11 pN. Finally, Lionnet et al. [2Bj examined 
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the unwinding kinetics of the gp41 hclicase. Like the T7 helicase, the unwinding velocity 
depends strongly on the value of the tension destabilizing the ss-dsDNA junction. From 
these experiments, we surmise that while different helicases exhibit very different unwinding 
velocity responses to force, the processivity always increases appreciably with force. Here, 
we provide a theory that explains the differing dependencies of unwinding velocity and 
processivity on force. 

Some aspects of the varied responses in hclicase velocities have been qualitatively justified 
in a previous work [3D], based upon an insightful theoretical model of helicases originally 
proposed by Betterton and Julicher (BJ) [DTI IDIfj . The model quantifies the crucial ideas 
of “active” and “passive” unwinding by a hclicase-a classification that has been the basis 
for understanding helicase mechanisms for a number of years [H El EH ESI EHl E1H3BJ . A 
passive hclicase utilizes the thermal breathing of the single strand-double strand junction of 
a nucleic acid, to opportunistically step in front. Since the double strand (ds) closes on itself 
on average, the hclicase frequently faces a barrier and hence moves ahead relatively slowly 
(as compared to its unimpeded single strand velocity). An active helicase on the other 
hand, promotes dsDNA melting by lowering the stability of base-pairs at the 
junction, either by utilizing free energy from ATP hydrolysis or the binding free 
energy of the protein at the junction. The resulting increase in unwinding velocity 
of an active hclicase depends on the extent of the destabilization of the ds junction due 
to interaction with the hclicase. This suggests that the helicase can be active to various 
degrees. 

An “optimally active” helicase would unwind at a maximum velocity that is close to or 
equals the single strand velocity, Vss [32] . For helicases with negligible back stepping rates, 
it was argued in [30] that when a hclicase is passive, an external force can assist in the 
opening of the double-strand junction, thus increasing the velocity of the hclicase. On the 
other hand, since the unwinding velocity of an optimally active hclicase is ~ Vss, an external 
force should not increase the unwinding velocity significantly. These physically motivated 
arguments suggest that the velocity of a passive hclicase should increase appreciably with 
force while the velocity of an optimally active hclicase will be similar to Vss- Manosas et 
al [30] did not investigate the processivities of helicases, and hence their work provided 
only an incomplete understanding of the nature of hclicase motion. By generalizing the BJ 
model [3TH33] to include force dependence for both velocities and processivities, we show 
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that even though the velocity can vary with external force in a hclicase-dependent manner, 
the unwinding processivity always increases significantly with force. This is a remarkable 
result, as it predicts a universal behavior — unlike the velocity, the processivity will increase 
rapidly with external force for all helicases regardless of their architecture. In addition, our 
work highlights another surprising hding — the more active a helicase is (the stronger its 
interaction with the single strand-double strand junction), the less processive is its motion. 
The effect of increasing the external force at a given interaction strength is 
therefore, opposite to the effect of increasing the interaction strength at a given 
external force. Finally, the sequence dependent behavior of the unwinding velocity and 
processivity show complex behavior depending on the percentage of GC content. We predict 
that details of the energy landscape of base pair opening, GC content, and the extent to 
which the helicase is active determine the unwinding velocity and processivity. 

II. MATERIALS AND METHODS 

As illustrated schematically in Fig. [lj in the BJ model the nucleic acid (NA) is represented 
by a one-dimensional lattice with n denoting the lattice position of the helicase, while m 
specifies the position of the ss-ds junction. Models similar in spirit, but differing sig¬ 
nificantly in formulation and application, have also been used to analyze motility 
of kinesin [39] . We extended the BJ model to include a constant external force or tension 
(F), which is applied to the complementary termini of the NA, in order to investigate the 
dependence of velocity and processivity on load. In accord with the single molecule optical 
tweezers experiments [28j 00], tension is applied in a manner that increases the opening 
rate of the junction and decreases the closing rate while maintaining detailed balance. 


If the helicase and ss-ds junction are in proximity, there is an effective interaction. The 
passive unwinding mechanism is realized by a hard-wall coupling potential: 



( 1 ) 


where j = m — n (Fig. [l|. Similarly, we represent an active unwinding mechanism by a 
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FIG. 1: A schematic illustration of the extension of the Betterton and Jiilicher model for helicases. 
The position of the helicase (black circle) on an underlying ID lattice representing the nucleic acid 
substrate is denoted by the variable n, while the variable m refers to the location of the ss-dsNA 
junction. At infinite separation between the helicase and ss-ds NA junction, n —> n +1 transitions 
occur at rate k + , while n —> n — 1 transitions occur at rate k ~. Similarly, m —> m + 1 transitions 
occur at rate a, and m —> m — 1 transitions occur at rate (5. The tension F is applied to the ends 
of the nucleic acid. 


coupling potential consisting of both a step-function and a hard-wall: 

°o (j < -1) 


U(j) = 


U 0 (j = 0) 

0 (j > 0 ), 


( 2 ) 


where the energy Uq is in units of ksT. 

Unwinding velocity of a helicase. The mean velocity (V) of the helicase/junction complex 
is proportional to a sum over all j of a product of the probability of being at separation j 
and the net rate at which the centroid coordinate (l = m + n, Fig. [I]) increases: 

v=l Yi di +- k - ft) ft, (3) 

j 

where ccj is the rate constant associated with junction opening when the helicase and junction 


5 











are at separation j. Similarly, fy, kf, and kj , are the j-dependent rate constants associated 
with junction closing, helicase forward stepping, and helicase backward stepping respectively. 

Since experiments are performed at a constant temperature and F, opening and closing 
of the ss-ds junction can only be attributed to thermal fluctuations or the applied tension. 
Thus, the ratio of the rate at which the junction opens (cty_i) to that at which it closes {fij) 
satisfies: 


a J - 1 = ® e [U(j-l)-UU)] e AG F 


Pi 


P 


a 

P = 6 


-AG 


(4) 


where a (/ 3) is the j —> oo junction opening (closing) rate, AG is the stability of the junction 
base pair in the absence of either force or helicase, U(j ) is the value of the coupling potential 
when the helicase and junction are at separation j , and A Gp is the destabilization in the 
free energy of the basepair at the junction, caused by the applied force F. For simplicity, we 
choose A Gp = FAx (see further discussions in the Discussion section), where Ax is roughly 
twice the length of ssDNA separating two bases. When ATP hydrolysis is tightly coupled 
to helicase transitions, the hopping rates approximately satisfy a relation akin to detailed 
balance I3T1 : 

h+ U+ 


k 


'j -1 


^p-iuu-v-um 

k- 


where k + = lim k+ and k = lim k j _ 1 . Note that when the helicase traverses on a single 
strand, there is no single strand-double strand junction confronting the helicase, and hence 
the hopping rates that describe its motion are k + and k~. Thus, the single strand velocity is 
Vss = k + —k~. For a number of helicases like T7, T4, NS3, UvrD which translocate 
directionally on ss nucleic acids, the forward rate k + can be assumed to be much 
larger than k~, near the physiological ATP concentrations (~ 1 mM, close to 
saturating conditions for most helicases) that are usually used in experiments 
[ ZTL 128H301 3BJ. We therefore work under this approximation throughout this 
article. While unwinding the double-strand however, the forward rate decreases 
and the backward rate increases (as evident from Eq. [5] and the equations for 
the individual rates below), making back-stepping much more likely. A high 
probability of backstepping was indeed observed in a recent experiment on the 
XPD helicase |37]. Finally, individual rates are affected by the coupling potential and 
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the applied tension as follows: k 7 = k + e~^ u ^~^~ u ^\ kj_ x = k-e^f-WU-^-uti)] ^ f3j = 
f3e~f [ C7 0' _1 ) _C7 h’)+ FAa; ], and aq_i = where / (0 < / < 1) is the 

location (in units of lattice spacing) of the transition state separating the closed and open 
states. Note that the applied tension F, affects only the nucleic acid breathing rates (acj-i 
and j3j), but not the hopping rates. 

The rates of a passive hclicase are independent of j for j > 1, and at j — 1 we have 
kf = (3i = 0. The force-velocity relation for a ‘hard-wall’ (passive) hclicase, calculated 
using Eq. [3j is: 


Vhw — 


a'k+ - p'k~ 
P' + k+ ’ 


( 6 ) 


where Vhw = Vhw(A,/), a' = ae~^~^ F/ ^ x and (3' = (3e~f FAx . For convenience, while 
referring to Vhw hi the rest of the text, we explicitly show the functional dependence of 
only the parameter pertinent to the particular discussion. Note that for F — 0, Eq. © 
coincides with Eq. (22) in reference [32]. By making the approximation k~ « 0, we obtain 
Vhw ~ j ( i + fc+/^ ) and vkwT^=o) ~ ( m+/? )- Fin ally, making the approximation 

(3' k + , we conclude that: 

V nw( F ) _ fax (7) 

Vhw (F = 0) ~ J 

(The approximation f3' k + is typically valid because f3 — 10 5 -10 8 s _1 , while k + ~ 1- 
10 3 s _1 ). A passive hclicase must wait for the junction to open to step forward and prevent 
the newly separated base pair from reannealing. From Eq. [7] it follows that the application 
of force exponentially increases the probability that the junction is open relative to the 
probability that it is closed, resulting in an exponential increase in the unwinding velocity 
relative to the F — 0 value. 

For a single step active hclicase (modeled with Eq. [2j) the rates are independent of j for 
all j > 1. A straightforward calculation of Eq. [3] leads to 


Vi = c'+ (l-d) e~f u0 

Vhw d + (l- d) e~ u ° ’ 1 j 

where Vj = V\{F, /, f/ 0 ), d = ( y ae~V~ 1 ') FAx _)_ k~)/{k + + f3e~V Ax ) and the subscripts 1 
and HW denote the step (active) and hard-wall (passive) coupling potentials. When the 
step potential goes to zero (f/ 0 = 0), the unwinding velocity is just equal to the hard 
wall (passive) hclicase velocity as can be seen from Eq. [8j As with Vhw, while referring 
to Vi in the rest of the text, we show only the functional dependence of the parameters 
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pertinent to the particular discussion. Finally, using the results of Eqs. [7] and [ 8 ] we find: 

W]B)) = ( L 7+ll-J)e~^ ) eFAX ( cHl-c]e : ^ ’ Wll6re C = C '( F = 0) = (« + k~)/(k + + /3). 

Helicase processivity. The processivity of a hclicase has been defined in various ways 
[ 33| : (1) the mean attachment time of the helicase (r), (2) the average number of base 
pairs unwound in a single binding event (5m) (in other words, the average number of base 
pairs by which the junction moves ahead before the helicase detaches), and (3) the average 
number of base pairs translocated before the helicase detaches (<5n); (5m) and (5n) (termed 
unwinding and translocation processivities respectively) in principle conld be very different. 
For instance, if the helicase binds very far from the junction (j 0 1), the double strand 

will close rapidly before the helicase can translocate by a significant amount, thus making 
the unwinding processivity ((5m)) negative. However, if the helicase binds very close to the 
junction (j 0 ~ 1 ), (5m) and (Sri) are almost identical [33]. The double strand would always 
have a larger closing rate as compared to the opening rate, thus making the strands close 
on average. This would hold even in the presence of external forces, as long as they are less 
than the force needed to unzip the double strand. In such physically relevant situations, 
after a very brief transition period the hclicase is likely to be very close to the junction, and 
hence the relevant initial condition can always be taken as jo — 1. 

In the single molecule experiments measuring unwinding velocity and processivity, typi¬ 
cally a DNA hairpin serves as a model double strand. The upstream single strand overhangs 
allow the hclicase to load |2T5. SOI G2E 3DJ. At forces less than the critical force required to 
unzip the hairpin in the absence of a hclicase, the arrival of a helicase at the junction and 
subsequent unwinding causes the end-to-end distance to increase. As a result, the presence 
of the helicase on the hairpin can only be discerned by observing the sudden change in the 
end-to-end distance, which happens when jo ~ 1. In this work therefore, we work with 
the initial condition j 0 = 1, and hence the unwinding processivity (5m) and translocation 
processivity (5n) are almost identical. 

To model a helicase with finite processivity, we incorporate an unbinding rate 7 j which 
depends implicitly on the separation j through the relation 7 j = [33]. We assume 

that U(j) —)■ 0 as j —» 00 . It is physically reasonable that the unbinding rate should 
increase as U(j) increases, because the repulsion between the helicase and the 
nucleic acid would cause the former to dissociate. This has indeed been observed 
in experiments — the dissociation rate of T7 from dsDNA is roughly 100 times 





larger | |21| than that from ssDNA [41 ], at the same dTTP concentration of 2mM. 
The exponential dependence of 7 ^ on U (j) is further justified by noting that it is 
built on the celebrated Bell model, which has been used successfully to describe 
unbinding of cell-adhesion moleclues and other complexes. Since the velocity of 
unwinding is unaffected by introduction of the unbinding rate 7 the analytical results 
derived above for the unwinding velocity hold good even for finitely processive helicases. 
For the physically relevant initial condition j 0 ~ 1 (explained above), the velocity of active 
unwinding is given by: 

Vj ~ (5m) / (r) ~ (5ti)/(t). (9) 

For the initial condition jo ~ 1, the two expressions for unwinding velocity given by Eq. [ 8 ] 
and Eq. [9] are equivalent. Betterton and Jiilicher [33] derived the following expressions for 
the three measures of processivity: 


<t> 

(5m) 

(5n) 


(SI) + (Sj) 
2 

(Si) - (Sj) 
2 


( 10 ) 


where the parameters {Rj} are obtained by solving the following infinite set of second order 
recurrence relations: 


- Sjjo ~ ~ + kj + ocj + fa + 7 j) Rj + (otj- 1 + k^) Rj- 1 + (f3j+i + kj~ + 1 ) Rj+i, ( 11 ) 


jo is the value of j at time t — 0, 5jj 0 is the Kronecker delta, (Sj) = .( j — jo)^jRj, and 


(51) is given by the expression: (51) = 6 (1 + a — y~ — ^ 


-1 r 


y- 


y+ 


(i-y -) 2 (i-y +) 2 


, where 


y± are the roots of the equation y 2 — (1 + a)y + (a — b) = 0, a = (1 + p)/q, b = 1/q 
and p = £j(aj + kJ~)Rj, q = J2j(Pj + k~)Rj. BJ derived Eq. |TTj by taking the Laplace 
Transform of the time evolution equation of P(j, l, t), the joint probability density of finding 
the helicase-junction system in state (j, l) at time t. P(j,l,t ) obeys a Master equation 
accounting for the movement of the helicase in the forward and backward direction as well 
as detachment from the NA. 

When the helicase is passive, U(j ) = 0 for j > 0 (Eq. [Tj) and hence the dissociation rate 
7 j = 7 for all j. As a result, the mean attachment time for a passive helicase (t)hw is given 
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by 


for all values of the external force F. 


III. RESULTS 


(t)hw — ~ , 
7 


( 12 ) 


In order to determine the velocity and processivity, we solve the sparse linear system 


given by Eq. 11 for {Rj}. All our results have been numerically obtained for the 
case where the back-stepping rate of the helicase in the absence of dsDNA (k~) 
is small compared to the forward rate in the absence of dsDNA ( k + ). We use 


a grid of size M — 10, 000 in j to solve Eq. 11 A larger grid size than that used by BJ 
(M = 100 was used in Ref. [33J ) is necessitated by the external tension applied to the NA 
substrate. From a numerical standpoint, as FAx/ksT increases, {Rj} converges to zero 
increasingly slowly, thus requiring more and more terms to guarantee the convergence of the 
sums used to calculate (r), (Sj), and (SI). For F Ax/ksT > 1.95, a grid of size M = 10 6 in 


j proves to be insufficient to solve Eq. (11). It is not possible to explore forces greater than 


1.95 ksT/ Ax because such large forces result in the melting of the duplex. This is because 
our choice of a = 10 5 s _1 and (3 = 7 x 10 5 s -1 correspond to AG = 1.95 k B T. At these 
forces, (8m) is very large and (Sn) ~ (r)(fc + — k~). These results are confirmed in Fig. [2j 
which shows (in addition to the numerical results) simulation results at these larger forces 
(see below for additional information). 

Unwinding velocity. To understand how the unwinding velocity depends on the model 
parameters, we plot the F-dependent unwinding velocity V\ (F,U 0 ) for several values of U 0 
(Fig. §• This figure reveals several interesting features. (1) Fig. [2] shows that for a small 
value of / = 0.01 < 1 (/ is the transition state location along the reaction coordinate 
separating the closed and open states of a base pair), if the helicase is passive (Uq = 0) or 
weakly active (Uo < AG), the unwinding velocity is highly sensitive to the external force. 
As Uo increases, making the helicase increasingly more active, the unwinding velocity is 
less sensitive to increase in force. In the parameter range / << 1 and Uo/ksT » 1, the 
helicase is “optimally active”. (2) The parameters used in Fig. [2] correspond to a single¬ 
strand translocation velocity Vgs=0.99 bp/s. In the optimally active regime, the unwinding 
velocity is close to Vss for a large range of forces. This observation along with point (1) 
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FIG. 2: Unwinding velocity V\(F, XJq) as a function of the tension (F Ax/ksT), for various coupling 
potentials Uo/ksT = 0 (red circles), 1 (orange squares), 2 (yellow diamonds), 3 (green downward 
triangles), 4 (green upward triangles), 5 (cyan rectangles). The solid black lines are numerical 
results obtained using Eq. [9] and Eq. Each symbol represents an average of 1000 independent 
Kinetic Monte Carlo simulations. For forces exceeding F Ax/ksT = 1.9, the duplex melts and it 
is no longer possible to numerically solve the system given by Eq. but the robustness of our 
simulation algorithm allows us to explore this regime confirming that the duplex has melted. The 
parameters used were / = 0.01, a = 10 5 s , (3 = 7 x 10 5 s , k + = lbp/s, k~ = O.Olbp/s, and 

7 = 0.01s" 1 . 

qualitatively describe the velocity behavior of a number of hclicases [30] • For example, UvrD 
is optimally active while T7 is passive or at best, only weakly active. (3) When U 0 = 0, the 
results revert to those expected for a hard wall coupling potential. Thus, as illustrated in Fig. 
[lj the effect of tensile force at Uq = 0 is simply to increase the mean velocity exponentially 
relative to the zero force value (Eq. [7]). 
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FIG. 3: Ratio of the mean attachment time of an active to a passive helicase (t)/(t) hw as a 
function of F Ax/ksT, for Uo/ksT = 1 (black solid line), 2 (orange dashed line), 3 (green short 
and long dashed line), 4 (blue short dashed line), 5 (red dotted line). (T)/(r) HW increases with 
FAx/ksT and decreases with increasing Uo/ksT. Interestingly, the mean attachment time is 
unaffected by / (all curves with the same value of Uo/ksT can be superimposed). Parameters 
were a = 10 5 s -1 , f3 = 7 x 10 5 s -1 , k + = lbp/s, k~ = O.Olbp/s, 7 = 0.01s -1 , 


used to solve Eq. 


11 


jo = 1, and M = 10 4 


Processivity and lifetimes. For finitely processive helicases, the primary results of numer- 
are plotted in Figs. [3j [4j [5] and [ 6 j illustrating several points worthy of 


ically solving Eq. 11 
note: 

1. The variations in the average lifetime of the helicase (r) (normalized by the lifetime 
of the passive helicase) with F, f and U 0 (Eq. 10) are displayed in Fig. [3j Evidently, 
a passive helicase (Uo = 0 ) has the largest lifetime. As Uq increases, lifetimes rnono- 
tonically decrease at any given force value. This is not surprising since the lifetime 
is controlled by the detachment rate 7 j, which increases exponentially with increase 
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in Uq. As F increases, the opening rate of the double strand becomes larger, thus 
increasing the probability of the helicase to find an open adjacent base whenever it 
steps ahead. This results in fewer occasions where the helicase has to pay the extra 
Uq energy to plough ahead, resulting in larger lifetimes. At forces close to the double¬ 
strand rupture force, the helicase always finds a clear path ahead and hence rarely 
interacts with the junction, making the lifetime increasingly approach the passive he¬ 
licase lifetime. Finally, Fig. [3] also demonstrates that the mean attachment time of 
the helicase relative to that of a passive helicase depends on both Uq and F, but is 
insensitive to the parameter /. 




FAx/k B T FAx/k B T 




FAx/k B T FAx/k B T 


FIG. 4: Mean unwinding processivity of an active helicase (relative to that of a passive helicase) 
as a function of FAx/ksT , for Uo/ksT =1, 2, 3, 4, 5, 10, 15, and 20. For all values of Uq, 
the processivity always increases with increasing tension destabilizing the ss-dsNA junction and 
decreases with increasing step height. Unlike the mean attachment time (Fig. [3]), the unwinding 
processivity is highly sensitive to the kinetic parameter /, further confirming that the processivity 


is likely to exert a strong influence over the kinetics of unwinding. Parameters used to solve Eq. 11 
were a = 10 5 s -1 , /3 = 7 x 10 5 s -1 , k + = lbp/s, k~ = 0.01bp/s, 7 = 0.01s -1 , jo = 1, and M = 10 4 . 
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2. Fig. [4] shows the mean unwinding processivity of an active helicase relative to a passive 
hclicase. Since the translocation processivity is very similar, we do not show a separate 
figure and all results discussed for (Sm), also hold good for (Sn). It is evident from the 
figure that similar to the lifetime (Fig. [ 3 ]), the unwinding processivity is also maximum 
for a passive helicase, and monotonically decreases as the helicase becomes more active. 
This interesting result can be physically understood as follows: (Sm) depends on the 
unwinding velocity and mean lifetime as (Sm) = V\ (F, Uq) (r). At a given force, as Uq 
increases, V\ (F,U 0 ) initially increases (see Fig. |6^,b) but (r) decreases (see Fig|3j Fig. 
[6^,f). The rate at which these two quantities increase/decrease determines the trend 
for (Sm). Our results show that the rate of decrease of (r) is faster than the rate at 
which V\ (F, Uq) increases. This can be most clearly seen in Fig. [6] (b and f) where the 
velocity and mean processivities have been plotted as functions of Uq respectively for 
/ = 0.25. (r) decreases faster than Vi(F, Uq) increases, and hence the overall result is 
that (Sm) decreases as a function of Uq regardless of the GC content. (Fig. [bjl). 

3. The unwinding processivity depends strongly on the parameter /. For / = 0.01 and 
F = 0, (Sm )/ (Sm) nw > 0.8 even when the coupling potential U 0 = 20 k B T (pink 
dot-dashed line in Fig. [5]). For / = 0.05 and F — 0, the coupling potential Uq must 
exceed 5 k B T for (Sm)/(Sm) nw to decrease to less than 0.8 (red dotted curve in Fig. 
[4]). Thus, the effects of the barriers to translocation and NA breathing on (Sm) are 
the principal determinants of the ‘velocity’ of a processive hclicase (see also sequence 
effects below). 

4. The equation (Sm) = Vi (F, Uq) (t) gives key insights into the force dependence of 
processivity. As was discussed above, for a passive helicase with Uq = 0, the velocity 
Vrw increases rapidly as the force increases (Fig. [ 2 ]). However, the lifetime of a passive 
hclicase (t)hw does not change and is constant at (t)hw — I /7 (Eq. [l2|. This 
immediately means that (Sm) will increase rapidly with increase in force (solid black 
lines in Fig. [ 5 ]). In the other extreme limit of an optimally active hclicase (Uo/ksT 1 
and / <C 1), our discussions above show that V\ will not significantly change with 
force (Fig. [ 2 ]). However, unlike the passive situation, (r) increases rapidly with force 
(Fig. §■ This result follows again from the discussion in the previous points—as a 
hclicase becomes more active, the increased interaction with the junction reduces its 
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FIG. 5: Plots of the mean unwinding processivity of a helicase as a function of the applied tension, 
for Uo/ksT =0 (passive),1, 2, 3, 4, 5, 10, 15, and 20. The processivity always increases with 
increasing tension destabilizing the ss-dsNA junction and decreases with increasing step height. 
We used the same parameters as in Fig. [4j 

lifetime. An external force reduces the probability of this interaction which in turn 
results in an increase in (r), the lifetime of the helicase. As a result, (8m) will increase 
rapidly as the external force increases (Fig. [5^). For intermediate values of U 0 and 
higher values of /, the increase in (8m) with force is a result of contributions from 
both the velocity as well as the lifetime. Thus, the phenomenon of rapid increase in 
unwinding processivity as F increases, arises due to very different reasons, depending 
on how active the helicase is. These arguments lead to the surprising prediction that no 
matter how active (or passive) the helicase is, the processivity is sensitive to external 
force. Our prediction, which we believe is universal for unwinding helicases, is borne 
out in the few experiments that have analyzed the variation of unwinding processivity 
over a range of forces [2HI 3D] • This behavior is to be contrasted with the dependence 
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of unwinding velocity on F, which varies significantly when the hclicase is passive, but 
less so for active helicases. 


Sequence effects. The quantitative insights obtained for hclicase velocity (Fig. [2]) and 
processivity (Fig. [5]) as a function of force prompted us to use a similar model to explore the 
effects of NA sequence. Helicase unwinding and translocation can be modeled as a discrete- 
state continuous-time stochastic process [42] • Instead of adopting Eq. 11 to include sequence 
effects, we used the Kinetic Monte Carlo (KMC) method 


to simulate the model. We 
chose KMC in this case because grid sizes larger than M = 10 4 would be needed to solve 
Eq. [lTJ which contributes to numerical stability problems. The transitions in the KMC are 
stochastically implemented, and a given trajectory is generated till the helicase disassociates. 

From detailed balance it follows that j| = e -AG , where A G ~ 2 is the sequence averaged 
free-energy per base pair in the absence of the helicase. To incorporate sequence effects we 
let AG = g(Sm ), where 8m = m — m o, mo is the initial position of the ss-dsNA junction on 
the lattice and the sequence-dependent function g is constructed from the nearest-neighbor 
parameters provided in Table 3 of SantaLucia et al. [44]. For example, if the junction is a GC 
base pair and the downstream pair is AT we assign A G(8m) = 1.46 kcal/mol ~ 2.43 /cbT. 
Following BJ [32j, we assume (3 to be constant and then assign a = f3e~ AG . As the sequence 


composition is varied, the average free energy per base pair ranges from AG = 1.7 when 


of the sequence used to investigate 40% GC content was 


the GC content is 0% to AG = 2.95 if the GC content is 100%. We investigated the NA 
sequences organized in a block copolymer fashion. For example, the infinitely repeating unit 

5' GGGGAAAAAA 3' 

3' CCCCTTTTTT 5'' 
We investigated sequences with fractional GC content varying from 0 to 1 in increments of 

0.1. All simulations were performed with initial separation % = 1 and with the junction 

located at the first GC pair of the repeating unit of the sequence. 

N 

The mean attachment time is calculated as: (r) = 1/A^T;, where N is the num- 

2—1 

ber of simulations (N = 1000 for every data point we collected), and t, is the at¬ 
tachment time in simulation i. The mean translocation processivity is calculated using 
(8n) = 1 /N'ff ( n{ — n® ), where n ° is the initial position of the helicase is simulation 

i=l ' ' 

i and n{ is its final position. Finally, the mean unwinding processivity is calculated as 

N , \ 

(5m) = 1/A Yf [ m i — mm, where m- is the initial position of the ss-dsNA junction and 

i=l ' ' 
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m{ is its final position. 

The results of our simulations investigating sequence dependence are provided in Fig. [6j 
and we note the following points of interest: (1) Not surprisingly, processivity (both (5m) 
and (t)) and V± = (Sm) /(r) decrease as the fraction of GC content is increased. (2) The 
decrease in attachment time with increasing Uq is similar at both / = 0.01 and / = 0.25. 
The behavior of (5m) with increasing Uq at / = 0.25, however, differs substantially from its 
behavior at / = 0.01. At / = 0.01 (5m) is essentially constant with respect to Uq. This leads 
to an eventual saturation of (5m)/(r) at large Uq. At / = 0.25, on the other hand, sequence 
effects are more pronounced. For example, (dm) decreases substantially with Uq when the 
%GC = 0 (Fig. |6j black solid line), but negligibly when %GC = 1 (Fig. [6j grey dashed 
line). This behavior leads to a pronounced peak in V\ when %GC = 0, but a very small 
peak when the fractional GC content is 1. Indeed, at %GC = 100, we see that (5m) can 
actually become negative leading to a negative V\ for hclicases with too strong a coupling. 
The surprising finding that velocity and processivity can be negative can be 
explained as follows. If the base ahead of the helicase is very stable, with a large 
AG, and if the helicase back-stepping rate (during unwinding) is comparable in 
magnitude to the forward stepping rate, it is possible that on average the helicase 
steps backward before the base in front opens from a thermal fluctuations. In 
such a situation, the velocity and processivity would be negative. It remains to 
be ascertained if this is realized in practice. 


IV. DISCUSSION 

Effect of ssDNA elasticity: To simplify our analysis, we used the Bell model A Gp = FAx, 
for the effect of force on the destabilization of the junction base-pair. It might be more 
accurate to use AGf = 2ylog (■^sinh(FV)) [29, 05], because a freely-jointed-chain model 
has been shown to be appropriate to describe single strand DNA elasticity [IS]. As shown 
in Fig. [7^, e AGp for both the models are similar, and hence using the simpler Bell model 
does not make any qualitative difference in our results. 

Effect of back-stepping rate: We analyzed the model for the case where the back-stepping 
rate of the helicase k~ is much smaller than the forward stepping rate k + , in the absence 
of the double-strand junction. Most molecular motors fall in this regime, where k + k~. 
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FIG. 6: The dependence of (r), (5m), and (6m) /(t) on the step-height ( Uo/ksT ) for varying 
amounts of GC content for the sequence given in the text. Simulations were performed at Uq/ ksT = 
0,0.5,1, 2, 3,4, 5, 6, 7,8,9, and 10. Each data point corresponds to an average over 1000 independent 
Kinetic Monte Carlo simulations. All quantities decrease with increasing %GC. Interestingly, the 
mean attachment time (t) is again very insensitive to f. (5m) is, however, sensitive to /, leading 
to disparate behaviors for (5m) /(t). When / = 0.25, (5m)/ (r) shows a very distinct maximum 
when %GC=0.0 but a very weak maximum when %GC=1.0. When / = 0.01, (6m) /(t) shows 
saturating behavior with increasing Uq. Thus, sequence plays a crucial role in determining the 
kinetics of unwinding. Parameters used in the simulations were a = 10 5 s _1 , /3 = 7 x 10 5 s _1 , 


k + = lbp/s, k = O.Olbp/s, 7 = 0.01s x , jo = 1. 
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For example, the ratio k + /k~ was measured to be about 221 in Kinesin [TF] and is expected 
to be large for helicases as well [3HJ. As long as this holds good, which has experimental 
support, all the results obtained in this work will remain valid. 

Mechanism for increase in processivity with force is different for passive and active heli¬ 
cases: From the recent single-molecule experiments on three helicases [23', [28], ! 40], we have 
surmised that the dependence of unwinding processivity on F may be universal (Fig. [7}o 
and c ). Irrespective of whether the unwinding velocity increases rapidly with force or not, 
the processivity seems to increase with the external force. In this work, we have provided a 
theoretical explanation of this behavior. Based on the theory, we predict that the sensitivity 
of processivity to force should indeed be a universal feature of all helicases, whether they 
are active or passive. Our argument hinges on the observation that the processivity of a 
hclicase is very well approximated by the product of two quantities, the unwinding velocity 
and the attachment time of the hclicase: (5m) ~ V\ (r). The origin of the universal increase 
of (5m) with F is dramatically different for passive and active helicases. We have shown 
that when the hclicase is passive, V\(— Vhw) increases rapidly with force while (r) HW stays 
constant, independent of the force. In contrast, when the hclicase is optimally active, V\ 
hardly changes as a function of force while (r) increases rapidly as the force is ramped 
up. Thus, in both these extreme situations of hclicase activity, the processivity (5m) shows 
significant variation as a function of force. This leads to the prediction that irrespective of 
the nature of interaction of the hclicase with the ds junction, the processivity should always 
increase as the force is increased. 

We note en passant the relationship, (5m) ~ V\ (r), can be used to anticipate the 
effect of ATP concentration on (5m). If for a particular helicase (r) is independent 
of ATP concentration, it follows that processivity is determined by V\ . So we 
expect that decreasing ATP concentration should decrease V], and hence (5m). 
Interestingly, this was observed in an experiment on RecBCD sometime ago j4§] . 

Is there a general role for partner proteins in the replisome?: Our demonstration that 
hclicase processivity always increases with F regardless of the underlying architecture of the 
enzyme, is likely to be relevant in vivo as well. It is becoming increasingly evident that 
helicases function most efficiently only in the presence of the macro-molecular machinery 
with which it interacts in the cellular milieu msg. Significant increase in the unwinding 
velocity and processivity of a number of helicases has been reported p2M2], when coupled 
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FIG. 7: (a) Dependence of e AGp on F for two models: A Gf = FAx (blue curve), A Gf = 
2ylog (-^sinh(F7)) (red curve). The parameters Ax = 0.594 nm, L = 0.6 nm/nucleotide and 
l = 1.3 nm were chosen such that for both models, the critical force (force at which AG = A Gf) 
is 13.5 pN, a typical value for DNA hairpins. To be consistent with the rest of our analysis, AG 
was chosen to be 1.95 ksT. (b) and (c) Experimental data suggesting a universal behavior of the 
unwinding processivity as a function of force. Velocity (blue) and processivity (red) data on (b) 
the T7 helicase [2B] and (c) the NS3 helicase m- The data shows that the unwinding velocity of 
the two helicases can be strongly or weakly dependent on external force. The processivity clearly 
increases as F increases, for both the helicases. 
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to partner proteins like the polymerase or single-stranded binding proteins (SSB). This 
naturally suggests a possibility that by partially destabilizing the ss-dsDNA junction, the 
partner proteins could be mimicking external forces that are applied in single molecule 
experiments. 

A number of explanations have been proposed in the literature, to account for the effect of 
partner proteins on helicases m ■ For example, dsDNA destabilization, prevention of back- 
slippage of the hclicase, specific interactions with the helicase leading to longer helicase 
lifetimes on the nucleic acid track, are some of the commonly invoked methods by which 
partner proteins are conjectured to increase the efficiency of helicases. Though it is likely 
that some or all of these mechanisms together contribute in increasing a helicase’s speed 
and processivity, there is evidence to suggest that partial destabilzation of the ss-dsDNA 
junction may be essential. It has been known for a long time that by binding tightly to 
ssDNA, SSB can destabilize the bases of dsDNA [53j[5l], at least the first few bases near the 
ss-ds junction [55] . A similar double-strand destabilization effect caused by the polymerase 
was shown to be critical in explaining single-molecule data on the T4 hclicase [52]. To 
explain the sometimes order-of-magnitude increase in velocity/processivity of helicases, it 
therefore seems likely that an essential requirement is at least partial destabilization of a 
few base pairs at the ss-dsDNA junction, by partner proteins. 

If partner proteins indeed mimic external forces resulting in destabilization of the ss- 
dsDNA junction, we would predict that they would have distinctly different effects on T7 
(weakly active) and NS3 (active) helicases. The effects should be similar to the force re¬ 
sponse (Fig. [7] b,c) - both the velocity and processivity of T7 hclicase should be increased 
by partner proteins, while only the processivity should be enhanced of the NS3 helicase. Re¬ 
markably, this has indeed been observed in experiments. Rajagopal and Patel [51] showed 
that SSB increase only the processivity of NS3 while Stano et al [50] and previous works 
[55] showed that both velocity and processivity of T7 hclicase increases due to the presence 
of a polymerase. 

It is worth re-emphasizing that we are drawing a parallel between the effect 
of partner proteins on the ss-dsDNA junction, and the base-pair destabilization 
effect (reduction in AG) due to external mechanical forces. We are not predicting 
the effect of partner proteins on the individual rates whose ratio determine AG 
(see Eq. [4]). Therefore, destabilization of the junction by partner proteins could 
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happen either by increasing the dsDNA opening rate a or by decreasing the 
dsDNA closing rate (3. The latter mechanism was shown to hold in a particular 
study of the effect of single strand binding proteins on dsDNA [57J, and hence 
is another example supporting our claim. 

Our work therefore suggests a very general role for the partner proteins associated with 
helicases in the replisome. Even though these associated proteins may or may not increase 
the unwinding velocity of a hclicase depending on the specifics of the hclicase’s interaction 
with the dsDNA, they should universally increase the processivity of the hclicase. This 
seems to be borne out in experiments, and the generality of our statement can be tested in 
future experiments. In a related context, it was shown earlier that the interactions 
between separate domains of the same helicase could potentially lead to a higher 
velocity than either individual domain [58]. This mechanism could also play a 
role in understanding the reason for increase in velocity of different oligomeric 
states of helicases or helicases interacting with partner proteins. 


V. CONCLUSIONS 

There has been no unified explanation for the recent experimental observations that the 
unwinding processivity of helicases is always enhanced by force while the velocity shows 
a multitude of helicase-dependent responses. By extending the theory of Betterton and 
Jiilicher to include force effects, we have provided an understanding of these varied observa¬ 
tions. The velocity of unwinding depends crucially on the extent to which a helicase is active. 
Optimally active helicases show little or no change in velocity while the velocity of passive 
helicases is highly sensitive to external forces that destabilize the double strand. In stark 
contrast, we predict that the unwinding processivity of a helicase should always increase 
with force, irrespective of how active or passive it is. The reason for this universal behavior 
of the processivity however, depends on the nature of the hclicase. Our prediction is general, 
and seems to be borne out in structurally diverse helicases like T7, UvrD and NS3. Future 
experimental measurements of load-dependent processivity can test our assertion that asso¬ 
ciated partner proteins may have coevolved with helicases to increase the processivity, not 
the velocity of unwinding nucleic acids. Finally, our sequence dependent simulations show 
that the velocity can display non-monotonic behavior that is strongly dependent on GC 
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content as well as saturating behavior at very small /. Many of our predictions await future 
experiments. 
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